Intrinsic double-peak structure of the specific heat in low-dimensional quantum 

ferrimagnets 
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Motivated by recent magnetic measurements on AsCu3(P04)4 (A = Ca, Sr) and 
Cu(3-Clpy)2(N3)2 (3-Clpy = 3-Chloropyridine) , both of which behave like one-dimensional ferri- 
magnets, we extensively investigate the ferrimagnetic specific heat with particular emphasis on its 
double-peak structure. Developing a modified spin-wave theory, we reveal that ferromagnetic and 
antiferromagnetic dual features of ferrimagnets may potentially induce an extra low-temperature 
peak as well as a Schottky-type peak at mid temperatures in the specific heat. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.30.Ds, 75.40.Mg 



Recent progress on the theoretical understanding of 
low-dimensional (low-D) quantum ferrimagnets deserves 
special mention in their long history of research. A min- 
imum of the susceptibility (x)-temperature (T) prod- 
uct has been known as typical of 1-D ferrimagnets |Q. 
Although the T~ ^diverging \F at low temperatures is 
reminiscent of the ferromagnetic susceptibility, it turns 
into the high-temperature paramagnetic behavior show- 
ing the antiferromagnetic increase. Recently an explicit 
sum rule j|] for the ferrimagnetic susceptibility has been 
found: Spin- (S, s) ferrimagnetic chains behave similar to 
combinations of spin-(S' — s) ferromagnetic and spin-(2s) 
antiferromagnetic chains provided S — 2s. An epochal 
argument H on the ground-state magnetization curves 
of low-D quantum magnets stimulated broad interest 
in ferrimagnetic chains QJ^] and ladders || in a field. 
Spin-(iS, s) ferrimagnetic chains were found to exhibit 
2s magnetization plateaus without any bond alternation 
. The discovery of a metal-oxide Haldane-gap antifer- 
romagnet Y 2 BaM0 5 [§,§] and its rare-earth derivatives 
i? 2 BaNi05 |l0| , [rT[| directed our attention to 2-D mixed- 
spin magnets. Their magnetic double structure, that is, 
the coexistence of gapless and gapped excitations, was 
well interpreted from the point of view of coupled ferri- 
magnetic chains fll2| |. Nuclear-magnetic-resonance mea- 
surements 1 13 on an ordered bimetallic chain compound 



NiCu(C 7 H 6 N20 6 )(H 2 0) 3 -2H 2 revealed a unique field 
dependence of the relaxation rate, T^ 1 oc H~ x l 2 , which 
was found to be indirect observation of ferrimagnetic dis- 
persion relations |p^| . 

In this article, we discuss another hot topic on the 
ferrimagnetic specific heat (C). Intertwining double- 
chain ferrimagnets ^3Cu3(P0 4 ) 4 (A = Ca, Sr) Jl|] 
and a ferromagnetic-ferromagnetic-antiferromagnetic- 
antifcrromagnetic bond-alternating chain ferrimagnet 
Cu(3-Clpy) 2 (N 3 ) 2 (3-Clpy = 3-Chloropyridine) @, 
which are illustrated in Figs. |l](b) and|l|(c), respectively, 
were both reported to exhibit a unique temperature de- 
pendence of the specific heat: A minimum at low temper- 
atures and then a noticeable increase toward a Schottky- 
type maximum. Ferrimagnets generally possess a ground 
state of macroscopically degenerate multiplet and there- 



fore an applied field may induce a double-peaked spe- 
cific heat p7| . However, such an extrinsic mechanism 
should be distinguished from the intrinsic thermodynam- 
ics. Magnetic measurements [jlj|,|l8],[l9| were also per- 
formed for these materials in an attempt to evaluate the 
exchange couplings but no parametrization reported so 
far is so consistent as to fit the magnetization, suscep- 
tibility, and specific heat consistently. Thus motivated, 
we calculate the specific heat for a wide class of 1-D fer- 
rimagnets and reveal its intrinsic and generic features. 



(a) 
(b) 
(c) 



S s 



J J 



s s s s Jaf Jaf Jf Jf 



FIG. 1. Schematic representations of the bimetallic chain 
compounds (a), the trimeric intertwining double-chain com- 
pounds (b), and the tetrameric bond-alternating chain com- 
pounds (c), which are described by the Hamiltonians (la), 
(lb), and (lc), respectively, where smaller and larger bul- 
let symbols denote spins s = | and S > |, while solid and 
dashed segments mean antiferromagnetic and ferromagnetic 
exchange couplings between them, respectively. 

The simplest quantum ferrimagnet in one dimen- 
sion consists of two kinds of spins S and s al- 
ternating on a ring with antiferromagnetic exchange 
coupling between nearest neighbors and a series of 
such family compounds were extensively synthesized 
by Verdaguer, Kahn, and their coworkers jij. One 
of their works is bimetallic chains of general for- 
mula ACu(pbaOH)(H 2 0) 3 -nH 2 (A = Ni,Co,Fe,Mn; 
pbaOH = 2-hydroxy-l,3-propylenebis(oxamato)) p0| , 
which are illustrated in Fig. [l](a). We first consider these 
bimetallic chain compounds in order to verify the valid- 
ity of our method, a modified spin-wave theory, as well 
as to understand typical features of 1-D ferrimagnets, 
and then proceed to the above-mentioned trimeric cop- 
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per phosphates and tetrameric copper complex, whose 
Hamiltonians can be written as 

H = J^2(S n ■ s n + s„ ■ S n+ i) , (la) 

n 

H = [J\{s n ,l ■ S„,2 + «n,2 ' S n , 3 ) 
a 

+ ^2(Sn,2 ' + s n.3 ' Sn+l,a)] > (lb) 

W = [Jaf(s„,1 • S n ,2 + «n,2 ' S„ >3 ) 

-Jf(«71,3 ' «n,4 + «n,4 ' «n+l,l)] ' ( lc ) 

respectively. In the following, we take N as the number 
of unit cells and set the length of the unit cell to unity. 

Assuming the Neel-like order and introducing the 
bosonic operators for the spin deviation in each sublat- 
tice, we can expand the Hamiltonian (|l^) as 

H = E cl&ss +H + H 1 +O{S- 1 ), (2) 

where E c \ ass = —2SsJN is the classical ground-state en- 
ergy, Ho gives the free spin waves, and Hi describes two- 
body interactions between them, which are, respectively, 
the 0(S 2 ), OiS 1 ), and 0(S°) terms. It may be an idea to 
diagonalize Ho and Hi simultaneously. However, the re- 
sultant dispersion relations are gapped and thus misread 
the low-temperature ferromagnetic features [pl| inherent 
in 1-D ferrimagnets, such as the T 1//2 -vanishing C and 
the T _2 -diverging x- Hence we propose another treat- 
ment [^|| of the up-to-O(S ) bosonic Hamiltonian. Ho 
is diagonalized as 

Ho = Eo + J ("k a t a k + wj^ft) , (3) 

k 

where E = J^2k\. UJk — + s )] 1S tne 0(S 1 ) quantum 
correction to the ground-state energy, and a\ and 0' k 
are the creation operators of the ferromagnetic and an- 
tiferromagnetic spin waves of momentum k whose dis- 
persion relations are given by lo^ = u) k ± (S — s) with 
uj k = [(S - s) 2 + 45ssin 2 (fc/2)] 1 / 2 . Using the Wick the- 
orem, Hi is rewritten as 

Hi = Ei - JJ2 {5u-a\a k + 5u+^/3 k ) 

k 

+ Wresid , (4) 

where the O(S ) correction to the ground-state en- 
ergy and those to the dispersions are given by 
Ei = -2JN[r% + r| + + \fs/S)rir 2 ] and 

5lo± = 2{S + s)(ri/u; k )sm 2 {k/2) + {r 2 /VSl)[u k ± 
(S - s)} with A = (2N)- 1 J2 k [(S + s)/u k - 
1] and r 2 = -N- 1 Y Jk (v r S~s/uj k ) cos 2 '{k/2), while 
the irrelevant one-body terms H mc \ = —J{S — 
s) 2 (ri/VS^)E k [cos(k/2)/Lo k ](a k (3 k + a\p\) and the 
residual two-body interactions 7i ros id are both neglected 



so as to keep the ferromagnetic branch gapless. This pro- 
cedure may be recognized as the perturbational treat- 
ment of Hi to Ho- The resultant Hamiltonian is com- 
pactly represented as 

w - E % + J E a l a * + "ivl fa) > ( 5 ) 

k 

with oj^ — u> k — 8uj^ and E g — E c \ ass + Eq + Ei . 
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FIG. 2. Noninteracting (dashed lines) and interacting 
(solid lines) spin-wave calculations of the dispersion relations 
and the specific heat for the alternating-spin chains which 
are described by the Hamiltonian (la) and Fig. 1(a). Cor- 
responding quantum Monte Carlo calculations (x) are also 
shown for comparison. 

In Fig. ||, the dispersion relations up to the order 
0{S 1 ) and 0(S°), ui k and u5j^, are compared with the 
numerical findings obtained thr oug h a recently developed 
quantum Monte Carlo scheme p3(| . We find two distinct 
branches of spin- wave excitations. The gapless one is 
made of the elementary excitations reducing the ground- 
state magnetization and is thus of ferromagnetic aspect, 
while the gapped one, enhancing the ground-state mag- 
netization, is of antiferromagnetic aspect. The ferromag- 
netic spin waves indeed exhibit a quadratic dispersion 
at small momenta. The O(S ) quantum correction has 
a significant effect on the antiferromagnetic spin waves, 
whereas the ferromagnetic ones look almost free from in- 
teraction. 

The core idea Q of the so-called modified spin-wave 
theory can be summarized as reliably describing thermo- 
dynamics even in low dimensions by introducing a con- 
straint on the magnetization. We demonstrate its fer- 
rimagnetic version p5[ for the alternating-spin chains. 
Constraining the total magnetization to be zero, Taka- 
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hashi [|26| obtained an excellent description of the low- 
temperature thermodynamics of 1- and 2-D Heisenberg 
ferromagnets. His idea that the thermal spin deviation 
should be equal to the ground-state magnetization may 
be replaced by 

(S + s)N=(S + s)Y^(nj;+n+)/w k) (6) 

k 

for our ferrimagnets, where n k = X^n- n+ f^Pkin" , n + ) 
with Pk(n~ ,n + ) being the probability of n~ ferromag- 
netic and n + antiferromagnetic spin waves appearing in 
the fc-momentum state. Equation (JsJ) claims that the 
thermal fluctuation should cancel the Neel-state stag- 
gered magnetization instead of the uniform one, in re- 
sponse to the replacement of the ferromagnetic exchange 
coupling by the antiferromagnetic one. Minimizing the 
free energy 

k <r=± 

+kBTY, Pk(n-,n+)lnP k (n-,n + ), (7) 

k n~ ,n+ 

with respect to Pu at each k under the condition (||) as 
well as the trivial constraints J2 n - n + Pk(n~ ,n + ) = 1, we 
obtain the free and internal energies at thermal equilib- 
rium as F = E B + n(S - s)N - k^T £ fc £ ct=± + n%) 
and U = E s + J J2 k X) CT =± "fc^fc > wnere the optimum 

distributions are given by n k = {e^* _ 
with a Lagrange multiplier \i obtained through Eq. 
(g). The thus-obtained specific heat, together with its 
interaction-free version, is shown in Fig. ^. In the case 
of {S,s) — (1, |), we compare the present results with 
the quantum Monte Carlo estimates. Considering that 
the conventional antiferromagnetic spin- wave theory does 
not work at all for the 1-D thermodynamics, the present 
calculations surprisingly well describe the overall temper- 
ature dependence, including the y/T initial behavior and 
the Schottky-type peak. When we employ the interact- 
ing spin waves, the most significant improvement is the 
correction of the peak position. Although they some- 
what underestimate the peak height, the Schottky-type 
anomaly is then correctly located, which is essential to 
reproduce the overall structure. The specific heat may 
be double-peaked provided the antiferromagnetic gap is 
much larger than the ferromagnetic band width, which 
is expressed in terms of the interacting spin waves as 
Wt_ 0J k =iT ~ ^k=Q- This condition is satisfied for 
S ^S> 2s and the ferrimagnetic chain of (S, s) = (|, |) 
indeed exhibits a double-peaked specific heat. In the fol- 
lowing, we argue the other mechanism for a varied tem- 
perature dependence, which is of topological origin and 
may thus be valid for homometallic 1-D ferrimagnets. 
Considering the poor convergence of quantum Monte 
Carlo calculations at low temperatures, the present mod- 
ified spin-wave scheme is one of the most reliable and 
feasible approach. 
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FIG. 3. Interacting spin-wave calculations of the disper- 
sion relations for the trimeric chains which are described by 
the Hamiltonian (lb) and Fig. 1(b), where the ferromagnetic 
and antiferromagnetic excitations are distinguishable shown 
by dotted and dashed lines, respectively. 
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FIG. 4. Interacting spin-wave calculations of the specific 
heat for the trimeric chains which are described by the Hamil- 
tonian (lb) and Fig. 1(b), where individual contributions of 
the ferromagnetic and antiferromagnetic spin waves are also 
shown by dotted and dashed lines, respectively, the sum of 
which is equal to the total drawn by solid lines. 

The dispersion relations and the resultant specific heat 
of the trimeric intertwining chains are shown in Figs. 
H and ||, respectively. The system is analogous to the 
alternating-spin chain of (S, s) — (1, |) but displays an 
additional flat band which is gapped but of ferromagnetic 
aspect. With decreasing J2 in comparison with Ji, a sec- 
ond anomaly appears at low temperatures. The present 
tool advantageously enables us to observe each contribu- 
tion of the distinct excitation bands. We find that due 
to the existence of the mid band, the double-peak struc- 
ture is limited to rather small ratios J2/ J\ . In the previ- 
ous experiments [l5| l , the minimum related to the second 
anomaly was indeed observed, but the measured temper- 
atures (1.8 ~ 15 K) were not low enough to estimate the 
maximum of the bump. A significant increase and the fol- 
lowing Schottky-type main peak at higher temperatures 
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were not explicitly shown either, where the magnetic and 
lattice contributions should carefully be separated. Our 
calculations fully motivate further measurements |Q and 
make possible detailed analyses of them. 
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FIG. 5. Interacting spin-wave calculations of the disper- 
sion relations for the tetrameric chains which are described by 
the Hamiltonian (lc) and Fig. 1(c), where the ferromagnetic 
and antiferromagnetic excitations are distinguishably shown 
by dotted and dashed lines, respectively. 
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FIG. 6. Interacting spin-wave calculations of the specific 
heat for the tetrameric chains which are described by the 
Hamiltonian (lc) and Fig. 1(c), where individual contribu- 
tions of the ferromagnetic and antiferromagnetic spin waves 
are also shown by dotted and dashed lines, respectively, the 
sum of which is equal to the total drawn by solid lines. 

The tetrameric bond-alternating chains are discussed 
m Figs. |and|. We again find a double-peaked specific 
heat of topological origin. The lower two bands con- 
struct the low-temperature bump, while the upper two 
contribute to the main peak at mid temperatures. It 
is the double contribution that makes the second peak 
much more noticeable in the tetrameric chains than in 
the trimeric chains. The susceptibility measurements 
p9| , p7t were well interpreted by setting Jf/Jaf to 0.5. 
However, such parameters neither induce any detectable 
second peak of the specific heat within the Hamiltonian 



( |lq ) nor fit the recent observations detecting the 
second anomaly around 0.5 K. We have tried the op- 
posite parametrizations Jaf < Jf- However, the third 
band lying much closer to the lower two makes the low- 
temperature anomaly much more conspicuous and sup- 
press the Schottky-type peak, ending up with complete 
discrepancy with the observations. Further measure- 
ments and more extensive analyses, for instance, taking 
account of exchange anisotropy, are necessary for the to- 
tal understanding. 

The trimeric and tetrameric chain materials have been 
measured by separate groups under their respective moti- 
vations. Interestingly, however, the unique temperature 
dependences of their specific heat potentially originate 
in the same mechanism. We have revealed the intrinsic 
double-peak structure of the specific heat in low-D ferri- 
magnets of topological origin. From this point of view, 
homomctallic or organic ferrimagnets should further be 
synthesized and measured, where the double-peaked spe- 
cific heat may generically be observed. 

This work was supported by the Japanese Ministry of 
Education, Science, and Culture and by the Sumitomo 
Foundation. The numerical calculation was done using 
the facility of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo. 
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